Fragile Signs of Criticality in the Nuclear Multifragmentation 
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Deviations from an idealized equilibrium phase transition picture in nuclear multifragmentation 
is studied in terms of the entropic index. We investigate different heat-capacity features in the 
canonical quantum statistical model of nuclear multifragmentation generalized in the framework of 
Tsallis nonextensive thermostatistics. We find that the negative branch of heat capacity observed 
in quasi-peripheral Au-fAu collisions is caused primarily by the non-generic nonextensivity effects. 
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The nuclear multifragmentation process is studied in 
the energetic collisions of heavy ions (HI). In these colli- 
sions, strongly off-equilibrium transient system is formed 
which equilibrates at the later stage of the reaction due 
to dissipative processes. Perfectly equilibrated system, 
as assumed in most theoretical descriptions of the multi- 
fragmentation decay of the hot residue, is most probably 
never attained. This would not be a serious problem 
if the nuclear fragmentation process does not show any 
sign of the 'criticality' Indeed, the nonextensiv- 

ity of weakly off-equilibrium finite systems may qualita- 
tively modify both the picture of the two phase coexis- 
tence and signatures of the critical behavior in small sys- 
tems Q. On the other hand, these nonextensivity correc- 
tions to Boltzmann-Gibbs statistical mechanics (BGSM) 
have no measurable effects on standard signatures of the 
equilibrium such as the particle/fragment kinetic energy 
spectra or angular distributions ||^. Neither the caloric 
curve nor the negative heat-capacity branch measure- 
ments both put forward as an evidence for the nu- 
clear liquid-gas phase transition, can be interpreted un- 
ambiguously. The nonextensive effects due to the long- 
range interaction/non-Markovian memory effects or the 
multifractal phase boundary conditions which are 
crucial only in the 'critical region' and cannot be reli- 
ably estimated, not only constitute an integral part of 
the physics of HI collisions but also provide an essential 
limitation to the understanding of the multifragmenta- 
tion based on the BGSM. The practical solution to this 
problem would be to use those models of excited nuclear 
matter which describe physics of limiting two phases out- 
side of the 'critical region' and to parameterize the crit- 
ical region in simple terms. It is the aim of this work to 
illustrate this problematic in the thermodynamic (canon- 
ical) model of the fragmentation which is extended to 
include the nonextensive effects in the framework of the 
generalized thermostatistics Q . 

A reasonable starting point could be any 'realistic' 
thermodynamic fragmentation model (for the list of ex- 
amples see 1^,^). This choice offers several advantages, 
such as correct quantum statistics, correct definition of 



fragment sizes and the fragment binding energies. The 
Coulomb and surface effects can be tuned by analyz- 
ing the observable quantities far outside of the 'criti- 
cal region' and fragment excitations can be included, 
if necessary. Several models of this kind had an un- 
questionable success in describing reaction products and 
their properties from the regime of particle evaporation 
at low excitation energies to the explosion at about 5 
- 10 MeV/nucleon The new class of nonexten- 

sive thermodynamical models, can be formulated in the 
framework of the Tsallis generalized statistical mechanics 
(TGSM) g. TGSM is based on an alternative definition 
for the equilibrium entropy of a system whose jth micro- 
scopic state has probability pi : 

g„^fc ^~^'^' =k ^^P^~^^P' ^ k>0 (1) 

q-1 q-1 

where q is the entropic index and the normalization con- 
dition 

Y.P^ = ^ (2) 

i 

is used to get the second equality in (|l|). The limit 
q ~ 1 corresponds to the BGSM. It is easy to verified 
that such general properties as non-negativity, concavity 
and so on are preserved by this new entropy definition. 
The main difference between BGSM and TGSM is in the 
non-additivity of entropy in the TGSM. Indeed, for two 
independent subsystems A, B, i.e. such that the joint 
probability of A + B is factorized into : pa+b = PaPb, 
the global entropy in TGSM : 

Sq{A + B) = Sq{A) + Sq{B) + (1 - q)Sq{A)Sq{B)/k , 

is not equal to the sum of the subsystem entropies. 

The entropic index q sometimes can be related in ex- 
plicit way to other basic quantities such as internal energy 
or free energy of the system. Recently the Tsallis defini- 
tion of entropy (^ was reinterpreted in terms of incom- 
plete information theory 0. The condition (||) means 
that all the possible physical states are well-known and 
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counted. But for complicated systems we often in prac- 
tice do not know all interactions, cannot find the exact 
Hamiltonian, the exact solution of equation of motion 
and the exact values of physical quantities. Therefore, a 
part of information is lost and the normalization (|^) is 
violated because the set of the countable states becomes 
incomplete. By redefining the real probability in Eq.(|^) 
to effective one as pi — > p^, one can keep the Tsallis defini- 
tion of entropy (note the 'invariance' of the numerator in 
Eq.(|l]) with respect to chosen normalization condition). 
Here the difference qi = q — 1 is related to the extra 
entropy of the system due to neglected interactions. In 
this case the essential question is not the preferential use 
of microcanonical ensemble (for which the entropy is ex- 
tremized with keeping only the normalization restriction) 
or canonical one (when an additional restriction for the 
energy conservation in average is applied) but rather the 
point which interactions are included effectively by in- 
troducing a free parameter q. The situation discussed 
we meet in the multifragmentation where effects of ther- 
mal and chemical non-equilibrium, slow expansion of de- 
caying system as well as possible multifractal behavior 
of interphase boundary are out of consideration in the 
BGSM. 

Our further consideration is based on the canonical 
multifragmentation model The canonical ensemble 
method in TGSM was introduced in |^ . The main ingre- 
dient of the nonextensive canonical quantum statistical 
model of nuclear multifragmentation is the expression 
for the fragment partition function : 



(3) 



where a and z are the fragment mass number and the 
fragment charge number, respectively. The fragment par- 
tition probability equals : 



Pp-(a, z) = [Ljg{a, z)] '[1 + qi f3 epia, z)]~^^^^ 



(4) 



where ep{a,z) = /2M ^ U{a,z) and /3 = l/T. In 
the limit qi — > 0, Eq.(^ recovers the familiar expres- 
sion : pp{a, z) = exp {—(3 ep{a, z))/uji(a, z). The internal 
energy [/, which includes the fragment binding energy 
and the fragment excitation energy, the temperature- 
dependent surface energy, and the Coulomb interaction 
between fragments in the Wigner-Scitz approximation, 
is parameterized as in Q. In the dilute gas approxima- 
tion , the partition function of a whole system can be 
written as : 



{a,z)] 



a,z Nn{a,z)\ 



(5) 



where the sum runs over the ensemble TiA,z of different 
partitions of A and Z of the decaying system : {fi\ = 
{iVft(l,0),7V^i(l,l),...,7Vft(A,Z)} and Nn{a,z) is the 



number of fragments {a,z) in the partition {n}. In this 
approximation, the recurrence relation technique [p|jTl|] 
can be applied providing exact expression for Qq [A, Z) 

Given the partition function, the mean value of any 
quantity is [0 : 



(6) 



In order to ensure the proper normalization of 
q— averages (^), it is better to work with the generalized 
averages M : 



=<o>, /<1>, 



(7) 



These normalized mean values exhibit all convenient 
properties of the original mean values. Moreover, the 
TGSM can be reformulated in terms of ordinary linear 
mean values calculated for the renormalized entropic in- 
dex : q* = 1 + {q — l)/q. In particular, the total average 
energy and pressure of the system become : 

£q = ^< N{a,z)>q.AZ <e{a,z)>q. (8) 

a.z 

Pq^^<N{a,z)>q*AZ <p{a,z)>q. (9) 

a.z 

where < e{a, z) >q and < p(a, z) >q are given by : 



< e{a,z) >q = - — 



d /I - [ujq{a, z)Y 
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(10) 



/A It) /I- K(a,z)]-«i\ , , 

-pdvj[ ) 

and the average multiplicity of (a, 2)-fragments in the 
fragmentation of system {A, Z) is : 



,w N / ,Qq{A-a,Z-z) 

< N{a,z) >qAZ^ (^qia,z)- 



QM^Z) 



(12) 



The heat capacity at a constant volume (= 
d£q/dT \vf) is: 

Cv = X! ^(a^; a'^') >g* < si"-, >q. x 

1.^ a' ,z' 

X < e{a' ,z') >q' + N{a,z) >q*AZ x 



[< £^(a,z) >q. - < e(a,z) >2 ]} 



(13) 



where : 



< A(az;a z ) >q=< N{a, z)N{a ,z ) >qAZ — 

- < N{a, z) >qAz < N{a ,z') >qAz (14) 



and 
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< N{a,z)N{a ,z ) >qAZ= 

Qq{A — a ^ a , Z — z — z ) 



ujy{a, z)LLig{a ,z ) 



Qq{A-a,Z-z) 



(15) 



The heat capacity at a constant pressure Cp (= d{Eq + 
VqVj)/dT \-p^) can be calculated using the relation : 
Cp-Cv ^ TVf KT {dVq/dT \ vff , where kt stands for 
the isothermal compressibihty —{l/Vf)dVf/dVq \t) ■ 



1 



-(3Vf 



^ ^ < A(az;a z ) 



X < p{a, z) >q* <p{a ,z) >q. + (16) 

E,r/ N I f d < p(a, z) >„. 

<iV(a,z) >,.^z - ' , y 9 



dVf 



and 



^1 



^ ^ < A(az;a z ) 



X < p{a, z) >q. < e{a ,z ) >q. + (17) 
1 (d < p{a, z) >q* 



+ ^(O'^) >q'AZ 



dT 



\Vf 



One should stress that all these thermodynamical quanti- 
ties are calculated exactly^ without using the Monte Carlo 
technique. 

The upper part of Fig. 1 shows the temperature de- 
pendence of the pressure for various entropic indices q 
in systems with = 100, 200 and 300 nucleons and 
Zf) = OAAq protons. In the bottom part, the temper- 
ature dependence of the inverse thermal compressibility 
1/kt is shown. Zero of 1/kt {dVq/dVf |t= 0) corre- 
sponds to the pole of Cp and defines the boundary of 
the two phase coexistence region. For q > 1, there ex- 
ists a region of temperatures where 1/kt is negative and, 
hence, Cp becomes negative between the poles. In the 
BGSM limit, 1/kt has zeros for Aq = 100,200, whereas 
in heavier systems these zeros appear only for g > 1. 
An essential part of the pressure and, hence, of 1/kt 
is the Coulomb term. The inverse compressibility 1/kt 
never vanishes when the Coulomb term is neglected 
Since the Coulomb contribution to the pressure and the 
inverse thermal compressibility decreases in the Wigner- 

— 1/3 

Seitz approximation roughly as Aq , this particular sig- 
nature may not be seen in heavy systems in the BGSM 
limit. Existing data do not allow yet to pin down the 
Aq— dependence of the 'criticality' signatures. Neverthe- 
less, Fig. 1 demonstrates how fragile is the Boltzmann- 
Cibbs equilibrium 'critical' behavior. Small increase of q 
above the BGSM limit leads to an upwards shift of the 



critical temperature Tc which, for the same value of q, is 
higher in heavier systems. All these important changes 
take place in a narrow range of temperatures around Tc, 
beyond which the fragmenting system closely follows the 
BGSM limit. 
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FIG. 1. The dependence of the pressure (the upper part) 
and the inverse isothermal compressibility (the lower part) on 
the temperature T is plotted for system of different sizes and 
different entropic indices q : 1.0 (the solid line), 1.0005 (the 
dashed line), 1.001 (the dotted line). The freeze-out volume 
Vf corresponds to p/ = Ao/Vf = po/4. The calculated values 
of l/kr are multiplied by factors 10 and 3 for Aq =100 and 
200, respectively. 

Fig. 2 present the heat capacities Cy and Cp as 
a function of the average excitation energy : E* — 
£q{T, Vf) - £q{T ^0,Vo), where Vf is the freeze-out vol- 
ume, Vq = Ao/po and po is the equilibrium density at 
T = 0. Cv is a smooth positive function of the excitation 
energy for all values of q. The peak of Cv{E*), whose 
position is associated with the critical temperature T^, 
becomes more pronounced for higher q. Fig. 3 compares 
the heat capacity Cp vs E* /Aq for systems of different 
sizes : Aq = 200 and Aq 300. In the BGSM limit, 
the negative branch of Cp is seen only for Aq < 200. 
With increasing Aq, its position moves towards lower ex- 
citation energies. It should be noted that the critical 
density pc in the nonextensive fragmentation model |^ 
is relatively high. The 'global' critical point {Vc,Tc,Pc) 
for Aq — 100 corresponds to Pc/po = 0.547,0.783 and 
0.925 for g = 1, 1.0005 and q = 1.001, respectively. For 
Aq = 200, the global critical point exists only for q = 1 
whose value of Pc/po = 0.904 is close to that obtained 
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in statistical multifragmentation model using the recur- 
rence relation technique [^. For g > 1, the system is 
always found inside of the two phase coexistence region 
for all p < pq. This is consistent with the picture of 
nonextensivity driven first order phase transition. 
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FIG. 2. The specific heat at a constant volume Cv (the 
upper part) and at a constant pressure Cp (the lower part) 
are plotted vs the excitation energy per nucleon for various 
entropic indices q in the system with Ao = 100 and Zo — 40. 

The description of nuclear matter in terms of the Van 
der Waals fluid (see also jl^ ) yields much lower crit- 
ical densities {pc « 0.3po)- In this model, the boundary 
of the coexistence region has a bell-like shape and the 
line V = Vf crosses it in a single point. Consequently, 
the negative branch of heat capacity is not seen. In the 
nonextensive fragmentation model , the boundary of 
the coexistence region is skew with the top tilted towards 
smaller V what allows for two crossings with the line 
V — Vf and leads to the negative branch of Cp. 

In conclusion, the phase transition in the statistical 
nuclear multifragmentation models tends to disappear in 
heavy systems due to the weakening of the Coulomb con- 
tribution. This effect can be compensated by the nonex- 
tensive features of entropy due to either long-range cor- 
relations/memory effects or the fractality of the liquid- 
gas interphase, which both tend to strengthen signatures 
of the first order phase transition. The application of 
nonextensive canonical statistical fragmentation model 
for the understanding of experimental 'caloric' curve 
data Q and the negative heat capacity data in the 
'critical' region, consistently indicates deviation from the 
BGSM picture of the phase transition and q > 1.0005. 
This tiny variation of q, which cannot be detected either 



in the particle/fragment kinetic energy distributions or 
in the angular distributions, have strong measurable ef- 
fects on the event- by- event energy fluctuations of parti- 
cles/fragments in the region of phase coexistence. Hence, 
the mass-dependence of the criticality signatures is deter- 
mined by a subtle competition between Coulomb repul- 
sive interactions and nonextensive features of the entropy. 
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FIG. 3. The specific heat at a constant pressure Cp is plot- 
ted vs the excitation energy per nucleon for various entropic 
indices q in the systems with Ao — 200, Zo ~ 80 (the upper 
part) and Ao — 300, Zo = 120 (the lower part). 

For q > I, the negative branch of Cp is seen both in 
light and heavy systems. The range of excitation ener- 
gies corresponding to Cp < increases with increasing 
Aq. However, in heavy systems the negative branch of 
Cp appears uniquely for q > I. Both extension and local- 
ization of the negative branch of Cp in quasi-peripheral 
Au+Au collisions at 35A.MeV [Q, closely resemble re- 
sults of nonextensive fragmentation model for q ~ 1.0005 
and ^0 — 200 (see Fig. 3). This suggests that the ob- 
served effect is caused primarily by the nonextensive fea- 
tures of the entropy. The position of singularity of Cp at 
higher excitation energies increases sensitively both with 
the entropic index q and with the source size. Correct de- 
scription of its experimental value [p^ by the nonexten- 
sive fragmentation model assuming maximal possible size 
of quasi-projectile source in this experiment {Aq < 200) 
I p^ , means that we have determined lower limit for q. 

There are many sources of nonextensivity in mcso- 
scopic systems. Some of them , e.g. the formation of 
liquid-gas (fractal) interphase [ [l5[, have been pointed out 
in the microcanonical studies | [l5| , p6t . Most of these ef- 
fects are non-generic and, moreover, they are hard to 
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quantify. This is a principal obstacle in the meaning- 
ful characterization of nuclear multifragmentation data 
in the 'critical region' in the idealized picture of BGSM. 
On the other hand, the same value of entropic index 
{q ^ 1.0005) seems to be consistent with both the caloric 
curve and the negative heat capacity Q data, in 
spite of completely different kinematical conditions in 
these measurements. Moreover, the excitation energy of 
higher singularity of Cp seems to be the same both in 
quasi-peripheral Au-|-Au collisions at 35A.MeV |Q and 
in central Xe-|-Sn collisions at 32A.MeV and agrees 
with q ~ 1.0005. This surprising universality remains a 
puzzle at present. 
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